arXiv.Org 

Mathematics - Numerical Analysis 



Domain decomposition schemes for evolutionary 
equations of first order with not self-adjoint operators 

Petr N. Vabishchevich 



Submitted to arXiv.org January 12, 2011 

Abstract Domain decomposition methods are essential in solving applied 
problems on parallel computer systems. For boundary value problems for evo- 
lutionary equations the implicit schemes are in common use to solve prob- 
lems at a new time level employing iterative methods of domain decompo- 
sition. An alternative approach is based on constructing iteration-free meth- 
ods based on special schemes of splitting into subdomains. Such regionally- 
additive schemes are constructed using the general theory of additive operator- 
difference schemes. There are employed the analogues of classical schemes of 
alternating direction method, locally one-dimensional schemes, factorization 
methods, vector and regularized additive schemes. The main results were ob- 
tained here for time-dependent problems with self-adjoint elliptic operators of 
second order. 

The paper discusses the Cauchy problem for the first order evolutionary 
equations with a nonnegative not self-adjoint operator in a finite-dimensional 
Hilbert space. Based on the partition of unit, we have constructed the op- 
erators of decomposition which preserve nonnegativity for the individual op- 
erator terms of splitting. Unconditionally stable additive schemes of domain 
decomposition were constructed using the regularization principle for operator- 
difference schemes. Vector additive schemes were considered, too. The re- 
sults of our work are illustrated by a model problem for the two-dimensional 
parabolic equation. 
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1 Introduction 

Domain decomposition methods are widely used for the numerical solution 
of boundary value problems for partial differential equations on parallel com- 
puters. Stationary problems pUfT^I^I^ are the most extensively studied in 
the theory of domain decomposition methods. Numerical algorithms with and 
without overlapping of subdomains are used here in the synchronous (sequen- 
tial) and asynchronous (parallel) organization of computations. 

Domain decomposition methods for unsteady problems are based on two 
approaches [H]. In the first approach for the approximate solution of time- 
dependent problems we use the standard implicit approximations in time. 
After that domain decomposition methods are applied to solve the discrete 
problem at a new time level. For the optimal iterative methods of domain 
decomposition the number of iterations does not depend on the discretization 
steps in time and space [311]. The iteration- free domain decomposition algo- 
rithms are constructed for unsteady problems in the second approach. In some 
cases we can confine ourselves to only one iteration of the Schwarz alternat- 
ing method for solving boundary value problems for the parabolic equation of 
second order [5J[7]. Special schemes of splitting into subdomains (regionally- 
additive schemes pSll27j l are constructed, too. 

Construction and convergence investigation of the regionally- additive schemes 
is based on the general results of the theory of splitting schemes [TUIIT^IHI] . 
The most interesting for practice is the situation where the problem operator 
is decomposed into a sum of three or more noncommutative not self-adjoint 
operators. In the case of such a multi-component splitting the stable additive 
splitting schemes are constructed on the basis of the concept of summarized ap- 
proximation. Additively-averaged schemes of summarized approximation are 
interesting for using on parallel computers. In the class of splitting schemes of 
full approximation |19| we highlight the vector additive schemes based on the 
transition from the single original equation to a system of similar equations 
[n[2ll3T]. Additive regularized operator-difference schemes are constructed in 
the most simple way for multi-component splitting |18ll23j . where stability is 
achieved via perturbations of operators of the difference scheme. 

Peculiarities of domain decomposition schemes result from the selection of 
splitting operators. To construct the operators of decomposition for bound- 
ary value problems for partial differential equations, it is convenient to use 
the partition of unit for the computational domain [51l51 fTBl[^[^[2^[55] . In 
the domain decomposition method with overlapping a separate subdomain is 
associated with a function with values lying between zero and one. Domain de- 
composition methods for unsteady convection-diffusion problems are studied 
in works |17ll20ll30] . In the limiting case the width of subdomain overlapping is 
equal to the discretization step. In this case the regionally-additive schemes are 
interpreted as the decomposition without overlapping of subdomains but with 
appropriate boundary conditions of exchange. Domain decomposition meth- 
ods for unsteady boundary value problems are summarized in the books [141 
[T^ . More recent studies are presented in the work [32]. In this case we use 
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different constructions for the splitting operators and for operators of the grid 
problem at a new time level. 

In this paper we construct domain decomposition schemes for the first 
order evolutionary equations with a general nonnegative operator in a finite 
Hilbert space. Decomposition operators are constructed here separately for 
the self-adjoint and skew-symmetric parts using the partition of unit in the 
appropriate spaces. Two classes of unconditionally stable regionally-additive 
regularized schemes are proposed. Vector additive operator-difference schemes 
of domain decomposition are considered. The paper is organized as follows. In 
section 2 there is formulated the Cauchy problem for the evolutionary equa- 
tion of first order and the corresponding a priori estimate of stability is de- 
rived. Decomposition operators are constructed in Section 3. Problems with 
the self-adjoint operator and skew-symmetric one are considered separately. 
Unconditionally stable regularized additive schemes of domain decomposition 
are constructed in Section 4, with the additive and multiplicative perturbation 
of the operator of transition to the new time level. Vector splitting schemes are 
discussed in Section 5. In section 6 we consider a model boundary value prob- 
lem for the two-dimensional parabolic equation along with the results of using 
different domain decomposition schemes. The main results are summarized in 
Section 7. 



2 The Cauchy problem for the first order evolutionary equation 

Let H he a finite-dimensional real Hilbert space of grid functions with the 
scalar product and norm (•, •) || ■ ||, respectively. Let a constant (independent 
of time t) grid operator A is nonnegative in H: 

d . , d 

- ' dt dt ^ ' 

and E is the identity operator in H. We search the solution of the Cauchy 
problem 

^+Au = f{t), 0<t<T, (2) 

u(0) = (3) 

Problem (HJ-Q results from a finite-difference approximation in space for 
the approximate solution of boundary value problems for partial differential 
equations. Similar systems of ordinary differential equations arise in using the 
finite-element method as well as in applying the finite-volume approach. Let 
us obtain the standard a priori estimate for problem ([I])-®. 

Multiply equation ([2]) by u scalarly in _ff . In view of ([T]) we obtain inequality 

i|lHP <(/,"). (4) 

Taking into account 

if,u)<\\f\\\\ul 
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from Q we have 

|lHI< 11/11. 

In view of the Gronwall lemma we obtain the desired estimate 



\u\\<\\u^\\ + / ll/WIM^, (5) 





which expresses the stabihty of the solution with respect to the initial data 
and right-hand side. 

The emphasis of our work is on constructing approximations in time for 
equation Two- level schemes will be considered. Let t be a step of the 
uniform grid in time and let j/" — y{t"), t" = tit, n = 0, 1, ■■■,N, Nt = T. 
Equation ^ is approximated by the two-level scheme with weights 

ri+l _ n 

^ ^ + A(ay"+i + (l-a)y") = ^", n = 0,l,...,7V-l, (6) 

T 

where, for example, (/?" = /{at"^^ + (1 — cr)<"). It is supplemented by the 
initial condition 

= (7) 

Difference scheme ([5]), © has approximation error 0{t'^ + {a ~ 0.5)r). 
Grid analog of ([5]) is the estimate at the time level 

< ||y"||+T||^"||, n = 0,l,...,7V-l. (8) 
Let us prove the following statement. 

Theorem 1 Difference scheme {IJ), ^ is unconditionally stable at a > 
0.5, and estimate 0) holds for the difference solution. 

Proof We write ^ in the form 

yu+i ^ _^ _^ arA)- V", (9) 

where 

S={E + arA)-^ [E - [1 - a)TA) (10) 
is the operator of transition to the new time level. From ^ we have 

M'+'\\ = \\S\\m+r\\iE + aTA)-'^^. (11) 

For the last term in the right-hand side of (|lip in the class of operators 
(dl), under natural conditions cr > we have 

||(i? + arA)-V"ll < 

We show that if a > 0.5 then for the nonncgative operators A the following 
estimate holds 

II^II<1- (12) 
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In the Hilbert real space H inequality is equivalent [9j to fulfilment of the 
operator inequality 

SS* < E. 

In view of ([TU]) this inequality takes the form 

{E + aTA)-^{E - (1 - a)TA){E - (1 - (j)tA*){E + arA*)-^ < E. 

Multiplying this inequality on the left by {E + aTA)~^ and on the right by 
(E + <7tA*)~^, we obtain 

iE-{l- a)TA){E - (1 - a)TA*) < {E + aTA){E + arA*). 

It follows that 

t{A + A*) + (0-2 _ (1 _ a)^yAA* > 0. 

This inequality holds for the nonnegative operators A with a > 0.5. In view 
of we have from ((TT|) required estimate ([5]). 



3 Operators of decomposition 

To better understand the formal structure of the domain decomposition oper- 
ators, we give a typical example. We consider a model unsteady convection- 
diffusion problem with a constant (independent of time, but depending on 
the points of a computational domain) coefficients of diffusion and convection. 
Convective transport is written in the so-called (see, for example, [21]) sym- 
metric form. Let in a bounded domain f2 an unknown function u(x, t) satisfies 
the following equation 

a— 1 ^ 

-Ea^(fcWa^)-/(-'^)' -e^' o<t<T, (13) 

a — 1 ^ \ ^ / 

where fc(x) > k > 0, x G Supplement equation with the homogeneous 
Dirichlet boundary conditions 

u(x, t) = 0, X e 5r2, < i < T. (14) 

In addition, we define the initial condition 

m(x,0) = u"(x), XGX2. (15) 

Let us consider a set of functions w(x, t) satisfying boundary conditions 
(jl4p . We write the unsteady convection-diffusion problem in the form of differential- 
operator equation 

du 

— +Au = f(t), 0<t<T. t>0 (16) 
dt 
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We consider the Cauchy problem for evolutionary equation (1161) : 

u(0) = u". (17) 

Operators of diffusive and convective transport are treated separately, so that 
in (dH) 

A = C + V. (18) 

For the diffusion operator we set 

a—l ^ 

On set of functions (fT4)) in TL = C2{fi) diffusive transport operator V is self- 
adjoint and positive definite: 

V = V* > K.6S, S = S{f2) > 0, (19) 

where £ is the identity operator in H. 

Convective transport operator C is defined by the expression 

a— 1 ^ ^ 

For any Wa(x) the operator C is skew-symmetric in T-L: 

C = -C*. (20) 

Taking into account representation ^TE\\ . it follows from p^ . (1^0]) that > 
in 

Domain decomposition schemes will be constructed via the partition of 
unit for the computational domain fi. Let domain f2 consists of p separate 
subdomains 

12 = r?i U X?2 u ... U Hp. 

Individual subdomains can overlap. With separate subdomain i?^, a = 1, 2, 
we associate function r/a{x.), a = 1,2, ...,p such that 

where 

p 

5^77<,(x) = l, xe/?. (22) 

Q = l 

In view of (l?T|) . we obtain from the following representation 
p 

^ = ^„=Ca+2?„, a = l,2,...,p, (23) 
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where 



Similarly (gD]), we have 



2?a=I?:>0, C« = -C:, a = 1,2,..., p. (24) 

In view of in splitting the following property holds 

Ac>0, a = l,2,...,p, (25) 

and the self-adjoint part of operator A is splitted into the sum of nonnegative 
self-adjoint operators, whereas the skew-symmetric one - into the sum of skew- 
symmetric operators. 

It is convenient to represent diffusive transport operator V as follows 

v = g*g, g^k^/^giad, g* = -divk^/^, (26) 

with g : TL ^ H, where H = {C2{^)y is the corresponding Hilbert space of 
vector functions. From this structure, for "Da, a = 1, 2, ...,p we obtain 

V,, = g*7j^g, a = 1,2,..., p. (27) 

Similarly, Ca, a = 1,2, ...,p can be represented as 

Ca = ^{riaC + Cria), a = l,2,...,p. (28) 



Representations ((27| . ((28|) for operators of diffusive and convective transport 
demonstrate us clearly the structure of operators in individual subdomains 
in the splitting based on (|2ip . (|22p and allow to verify fulfilment of A 
similar consideration can be given for the operator of problem ©, ([3]). 

Let us divide the operator A into the self-adjoint and skew-symmetric parts: 

A = C + D, C=^{A-A*), D = ^{A + A*). (29) 

For the nonnegative operator D the following representation holds 

D = G*G, (30) 

where G . H H. For the partition of unit of the computational domain we 
consider the corresponding additive representation of unit operators E and E 
in spaces H and H , respectively. Let 

p 

Y^Xc^E, xa>0, a = l,2,...,p, (31) 

a=l 
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Y^Xc^E, xa>0, a = 1,2,..., p. (32) 

a=l 

By analogy with we use the sphtting 

p 

A=^A„, A^>0, a = l,2,...,p, (33) 

where 

A^^Co^ + D^, D^^DI>0, C^^-Cl, a = 1,2,..., p. (34) 

On the basis of ([5^ for the terms of the self-adjoint part of the operator A we 
set 

Ca = G*XaG, a = 1,2,..., p. (35) 
Decomposition of the skew-symmetric part is based on pip : 

Co.^\{XaC + Cxc.), a = 1,2,..., p. (36) 

Such an additive representation is a discrete analog of (l?fl) , (^5]) and is inter- 
preted as the corresponding version of the domain decomposition. 

4 Regularized domain decomposition schemes 

For the approximate solution of the Cauchy problem for equation ([2]), (O under 
condition (j33p we apply different splitting schemes. Transition to the new time 
level is based on solving p separate subproblems with individual operators 
Aa, a = 1,2,..., p. Taking into account the structure of the operators (see 
(p4)) -(l36 l) l we can say that these splitting schemes are regionally-additive and 
iteration-free. 

Currently, the principle of regularization of difference schemes is considered 
as the basic methodological principle for improving difference schemes [13) . 
The construction of unconditionally stable additive difference schemes [19] via 
the principle of regularization will be implemented in the following way. 

1. For the initial problem there is constructed some simple difference scheme 
(the producing difference scheme). This scheme does not satisfy the nec- 
essary properties. For example, in constructing additive schemes the pro- 
ducing scheme is not splitting one, or it is conditionally stable or even 
absolutely unstable. 

2. The difference scheme is written in the form for which the stability condi- 
tions are known. 

3. Quality of the scheme (its stability) is improved via perturbations of oper- 
ators of the difference scheme with preserving possibility of its computa- 
tional implementation as an additive scheme. 
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Concerning to problem © it is natural to choose as the producing 
scheme the following simple explicit scheme 

n+l _ n 

^ ^ + Ay" = ^", n = 0,l,...,7V-l, (37) 

T 

which is supplemented by initial conditions Q. Stability of scheme (|37)) is 
provided (see the proof of Theorem [T]) by fulfillment of inequality 

A + A*~ tAA* > 0. (38) 

Inequality psp with D > imposes appropriate restrictions on the time step, 
i.e. scheme ([25]) . ([57| is conditionally stable. Note also that if _D = than 
scheme ((29|) . ((37)) is absolutely unstable. Taking into account splitting ((33|) . 
we refer this scheme to the class of additive schemes. 

To construct additive schemes, we can take more general scheme (O, (O as 
a producing one. It is unconditionally stable at cr > 0.5. In this case the per- 
turbation of scheme operators is oriented only to receive the additive schemes 
preserving the property of unconditional stability. 

Regularization of a difference scheme in order to improve the stability re- 
striction (construction of a splitting scheme) can be performed via some per- 
turbation of the operator A. The second possibility is related to perturbation 
of the operator at the difference derivative in time (for our scheme ( ref 37) it 
is operator E). In constructing additive schemes, it is convenient to operate 
with the transition operator S, rewriting producing scheme (I57| as follows 

y"+i S'y"-|-r(p", n = 0,l,...,iV- 1. (39) 

For ((37)) we have 

S = E-tA. (40) 

The regularized scheme is based on perturbation of the operator S and has 
the following form 

= + n = 0,l,...,iV-l. (41) 

Let us consider general restrictions for S. 

To preserve the first order approximation, which has generating scheme 
((39)) . ((40)) . we subordinate the selection of S to the condition 

S = E -tA + 0{t^). (42) 

Stability of scheme ((^T)) in the sense that estimate ((5]) holds, is provided by 
the inequality 

ll^ll < 1- (43) 

In addition, the regularized scheme must be additive, i.e. the transition to 
the new time level is implemented via solving the individual subproblems for 
operators Aa, Q!=l,2,...,pin decomposition p3D . 
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The first class of regularized splitting schemes is based on the following 
additive representation for the transition operator of the producing scheme 

1 ^ 

S=-y^Sa, Sa^E-prAa, a = l,2,...,p. 
p ^ 

a — l 

The similar additive representation we also use for the transition operator of 
the regularized scheme 



S=-J2Sc., a =1,2,..., p. (44) 



a=l 



Individual terms Sa, a = 1,2, ...,p are constructed via perturbations of oper- 
ators Aa, a = 1,2, ...,p. By analogy with dTO)) we set 



S^ = {E + apTA,,)-^{E-{l-a)pTA^), a = 1,2,..., p. (45) 

If (7 > 0.5 (see proof of Theorem [T]) we have 

||5„||<1, a = l,2,...,p. 

In view of (HU it provides fulfilment of stability conditions P5)) . 
Using the representation 

Sa = E ~ pt{E + apTAa)^^Aa, a = 1,2, ...,p 

we can rewrite regularized additive scheme (j41[) , , (HSl) as follows 

n+l n P 

^ ~^ +^(g + apT.4^)~M^y"^y", n = 0,l,...,iV-l. (46) 

a— 1 

The comparison with producing scheme (I33p . (I37p shows that the regulariza- 
tion is provided by perturbation of the operator A. Our consideration results 
in the following statement. 

Theorem 2 Additive difference scheme^, Ii41\ l, |^5[ ) is unconditionally 

stable at a > 0.5, and for the numerical solution estimate of stability ^ with 
respect to the initial data and right-hand side holds. 

Numerical implementation of scheme (O, ([^S]) can be conducted as follows. 
Assume 



In this case we obtain 



Ji+l _ n 

^ ^— + {E + apTAa)-'Aay'' = ipl, a = 1,2,..., p. (47) 

pT 
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for the individual components of the approximate solution at the new time 
level ya~^^, a = 1,2, ...,p. Scheme (|T7)) can be rewritten as follows 

„,n+l _ „,n 

^ ^ + A^y-{<jy:+' + (1 - = {E + aprAM- 



In this form we can interpret scheme (j47p as a variant of the additively- 
averaged scheme of component- wise splitting [19 . 

The second class of regularized splitting schemes is based on using not 
additive (see (1441) ) but multiplicative representation of the transition operator: 

5=J|5o. a = 1,2,..., p. (48) 

Taking into account (HU, we have 
p 

5= []5„ + 0(r2), So,=E-tA^, a = 1,2,..., p. 

a=l 

Similarly ps)) . we set 

So, = {E + aTA^)-\E-{l-a)TA^), a = 1,2,..., p. (49) 

Under the standard restrictions a > 0.5 regularized scheme (|1T|) . (021) is 

stable. 

Theorem 3 Additive difference .scheme i41\ l, i48{ l, |/^ff[ ) is uncondition- 
ally stable at a > 0.5, and estimate of stability with respect to the initial 
data and right-hand side is valid for the difference solution. 

We present now some possible implementation of the constructed regular- 
ized scheme. Let us introduce auxiliary quantities 

yu+a/p^ a = 1,2,..., p and 
taking into account (|H|), (|48l) . we define them from the equations 

yu+a/p ^ a = 1,2, ...,p - 1, 

= 5py"+(P-l)/P + Tip''. (50) 

Similar to (|Tf)) we obtain from (15(11) 

y- ^ + {E + arA^r^A^y^+^'^'^yP = ^l, (51) 

r 

where 

;0, a = l,2,...,p-l. 



ip-, a = p. 



Rewrite scheme (1511) as follows 



..n+a/p _ n+{a — l)/p 

^ ^ + A„(a2/"+"/P + (1 - fT)y"+("-i)/f ) = ^l, (52) 
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where 

^l = {E + aTA^)^l, a = 1,2,..., p. 

Scheme ((52|) is a special variant of the standard component-wise splitting 
scheme [TOlfTSllBl] . But unhke these schemes of summarized approximation 
we constructed here the regularized schemes of full approximation. Regular- 
ized schemes (|1T|) . (H^ . (H5|) . constructed using additive representation pij) 
for the transition operator, are more suitable for parallel computations in com- 
pare with regularized schemes (j41L (|48)) . (l49t which are based on multiplicative 
representation (|48|) . 



5 Vector schemes of domain decomposition 

Difference schemes for unsteady problems can often be treated as appropriate 
iterative methods for the approximate solution of stationary problems. Great 
opportunities in this direction provide the vector additive schemes [11131]. 

Instead of the single unknown u{t) we consider p unknowns Ua, a = 
1, 2, which are determined from the system 

-^ + Y.Apup = f[t), a = l,2,...,p, Q<t<T. (53) 

p=i 

The system of equations (|53p is supplemented with the initial conditions 

u^{Q) = u\ a = l,2,...,p, (54) 

which follow from ([2|). Obviously, each function is a solution of problem ([2]), 
(|33p . The approximate solution of ©, ([3]), will be constructed on the 
basis of one or another difference scheme for vector problem ((55)) , ([51)1 . 
To solve problem (j53L (|M)) . we use the following two- level scheme 

— — + ^/3y;3 + 2^ ^fiVf) = V ' 

^=1 I3=a+1 

a=l,2,...,p, n = 0,1,...,A^- 1. (55) 

For this difference scheme we use the initial conditions 

2/a(0)=uO, a =1,2,..., p. (56) 

Numerical implementation of this scheme is based on the successive inversion 
of operators E -\- rAa, a = 1, 2, 

Theorem 4 Vector additive difference scheme \33l\) . ^56\) is uncondi- 

tionally stable, and for the components of the difference solution the following 
stability estimate with respect to the initial data and right-hand side 



k=l 



i^k _ 1 
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a = l,2,...,p, n = 0,l,...,N -1, 



(57) 



is valid. 

Proof To study vector scheme ([55)) . it is convenient to use the approach 
from the work [22 . Subtracting the a-th equation for y2_i_i from the a + 1-th 
equations of system ((55)l for y^+i^ S^^ 



{E -f tA^+i 



,,"+1 _ ,,n ,/"+! „ 

^ !Ja + l yg+l _ Ua Vg 

T T 

,,n+l 



a = l,2,...,p- 1. (58) 



Similarly, considering the equations for y" and , we obtain 



,,"+1 _ „n y" — _ (/^"-l 



(59) 



Taking into account nonnegativity of operators j4q,, a = l,2,...,p, from 
we obtain 



,,"+1 _ 



< 



a = l,2,...,_p- 1. 



(60) 



Similarly, from (j59p we have 



,,"+1 _ „« 



< 



(61) 



From ([50)1 . ([5T|) . we derive at each time level the following estimate 



,n „,n — 1 



< 



2/a - 



a = l,2,...,p, n = l,2,...,A^-l. 



(62) 



From (|62|) we get 



yS+' - yS 



< 



yi - 



a = l,2,...,p, n = 1,2, A^- 1. 



(63) 

From (|55p with a = 1, taking into account splitting (I33p and initial conditions 
(ISSl)!. we obtain 



yi-y? 



In view of (|S(7| we can rewrite inequality as follows 



y^' - yS 



<||(p°-A^.O||+^r 



/c=l 
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a = l,2,...,p, n = 1,2, A^- 1. 
Taking into account the obvious inequality 



(64) 



i2/r'ii< 



a = l,2,...,_p, 



we obtain from ([M]) required estimate ([57]) . 



We emphasize that the above stabihty estimates ([57| are received for each 
individual component y^^^, a — 1,2, Each of them or their linear com- 
bination 



^„+i ^ ^ Co,yl+^, c„ = const > 0, a = 1, 2, 



a=l 



can be treated as an approximate solution of our problem © , , ([55]) at time 
moment t — 



6 Model problem 

To illustrate possibilities of the constructed here domain decomposition schemes, 
let us consider the simplest boundary value problem for the parabolic equation. 
We consider the problem in a rectangle 

= { X I X = {xi,X2), < Xa <lai Q; = 1, 2}. 

In f2 the following boundary value problem 



Q — 1 " 

u(x, 0, X e < < < T, (66) 

u(x,0) = u°(x), xe/? (67) 

is solved. 

The approximate solution is searched at the nodes of a uniform rectangular 
grid in f2: 

w = {x I X = (a;i,a;2), x^^iaha, ia ^ 0,1, N^, Nah^ ^ la} 

and let uj be the set of internal nodes (uj — ujU du)). For grid functions ?;(x) = 
0, X G duj we define the Hilbert space H — ^2(1^) with the scalar product and 
norm 

{y, w) = Y^ y(x)w(x)/ii/i2, ||y|| = {y, yf^. 
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Approximating problem (j65p . (|66p in space, we obtain the differential- 
difference equation 

^ + Ay = 0, X € cj, < i < T, (68) 

where ^ 

Ay = -j^ivixi + hi,X2) - 2y{xi,X2) - y{xi - /ii, ^2)) 

- ^iyixi^X2 + - 2y{xi,X2) - y{xi,X2 - 112)), X e w. (69) 
In the space H the operator A is self-adjoint and positive definite [T?lll5] : 

A^A*>i6i+62)E, ,5„ = ±sin2^, a = 1,2. (70) 

Taking into account (|67)) . we supplement equation ((69|) with the initial condi- 
tion 

y(x,0) =M°(x), xew. (71) 

For simplicity, operators of domain decomposition in the investigated prob- 
lem ([S5|) - ([7T|) will be constructed without the explicit definition of operators 
G and G as well as space H, focusing on decomposition (I?!]) . 1^7^. We set 

^<^y = "/^^'^"(^i + 0-5hi,X2){y{xi + hi,X2) - y{xi,X2)) 



1 

hi 
1 
hi 



+ 72-'7a(^i ^ 0.5/ii,X2)(y(a;i,a;2) - yixi - hi,X2)) 



-^ria{xi,X2 + 0.5h2)iyixi,X2 + /12) - yixi,X2)) 



+ j^r]aixi,X2 - 0.5h2){y{xi,X2) - y{xi,X2 - h2)), a = 1,2,..., p. (72) 

In view of (HI]), dm we have 

p 

A = ^^,, A^^Al, a = 1,2,..., p. (73) 



Thus, we consider the class of additive schemes ([33)) . 

Numerical calculations for problem (j65l) - ((67|) are performed in the unit 
square {h = I2 = 1) where the solution has the form 

m(x, t) = sin(ni7ra;i) sin(n27rx2) exp(— 7r^(n^ -|- n2)t) (74) 

for a natural ni and 7^2. For this solution we set the corresponding initial con- 
ditions (|67p . Decomposition is performed with respect to one of two variables 
into four subdoniains (see Fig. [1} with overlapping. Disconnected subdomains 
can be considered as some single subdomain and the decomposition in Fig. [T] 
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Fig. 1 Domain decomposition 

can be treated as the decomposition into two subdomains described via two 
functions: rja = ria{xi), a = 1,2. 

For problems of type (|65|) - (l67t two cases of domain decomposition are high- 
lighted: decomposition with and without overlapping of subdomains. Methods 
without overlapping of subdomains are connected with the explicit formulation 
of certain conditions at the common boundaries. In our case a special problem 
at the interfaces is not formulated, but for algorithms without overlapping we 
can derive the corresponding exchange boundary conditions. 

For the domain decomposition methods the fundamental issue is exchange 
of calculation data between different subdomains. Standard explicit schemes 
can be used. In this case the domain decomposition can be associated with 
separate subsets of grid nodes:ajQ,, a = 1, 2, where uj — ujiU uj2- In the case of 
(|55|) - ([S7|) (the seven point stencil in space), the transition to the new time level 
via the explicit scheme for finding the approximate solution on grid Wq, a = 
1, 2 is performed using the solution values at nodes adjacent to the interface. 
We need to transfer the data of size ~ dwa, a — 1,2. For the approximate 
solution of problem (j68P "(|7i p we can consider two possibilities with minimum 
overlapping of subdomains. The first employs the domain decomposition with 
interfaces at integer nodes — the boundary nodes belong to several subdomains 
(two in our case of decomposition with respect to one variable). The second 
possibility is realized when the boundary of subdomains passes through the 
half-integer nodes of the corresponding variable. 

The variant of domain decomposition with boundaries through integer 
nodes is shown in Fig. [51 Assume that the decomposition is carried out in 
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!7i Q2 
Vi{9) mid) 



Fig. 2 Decomposition through integer nodes 



spatial variable xi, i.e. 6 = xi. The boundary of subdomains here passes 
through the node — 6i. Thus, for this decomposition operators ([7^ take the 
form ^ 

^^V = T2^y{xi,X2) ~ y{xi - hi,X2)) 



--^{y{xi,X2 + h2) - 2y{xi,X2) - y{xi,X2 - /12)), 



1 

^2 



^22/ = -r2(y(-^l + fT'li^-i) - yixi,X2)) 



1 

--^{y{xi,X2 + ^2) - 2y{xi,X2) - y{xi,X2 - /12)), xi = 9i. 

This decomposition can be associated with using the Neumann boundary con- 
ditions as the exchange boundary conditions. The relationship between the 
individual subdomains is minimal and requires exchange of data aX 6 ~ 6i. 
This case can be identified by the operators of decomposition as follows: 

i?(Xa) = [0,1], a = 1,2,..., p. (75) 

The values of rja{xi ± 0.5/ii, X2), ria{xi,X2 ± 0.5/ii), a = 1,2 for (17^ . ([71]) is 
equal to or 1. 
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Fig. 3 Decomposition through half-integer nodes 

The second possibility, which is associated with decomposition through 
half-integer nodes, is depicted in Fig. [31 In this case instead of (|75p we have 

i?(x„) = [0,1/2,1], a =1,2,..., p. (76) 
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At node 6 — 9i we use the difference approximation with the flux reduced by 
half. For the decomposition in the variable xi the operators of decomposition 
d72t seem like this 



Aiy = -^{y{xi,x2) - y{xi - hi,X2)) 



-7T2(y(^i'^2 + ^2) - 2y{xi,X2) - y{xi,X2 - /12)), 
^2?/ = --^{yixi + hi,X2) - y{xi,X2)) + ^(y(a;i,a;2) " yixi - hi,X2)) 



3 

'4^ 



(y(a;i,a;2 + ^-2) - 22;(xi,X2) - y{xi,X2 - /12)), a; 



'2 

For calculations in the domain !7i (see Fig. [3])) we employ adjacent to the 
interface data from the domain Q2 — at node theta = thetai. Thus, for this 
domain decomposition exchanges are minimal and coincide with the exchanges 
in the explicit scheme. 
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9i~l 9i Si+l/2 di+i 

Fig. 4 Decomposition through integer nodes with the width of overlapping 3h 



The considered variants of decomposition ([75)1 . ([75]) correspond to the mini- 
mum overlapping of subdomains. At the discrete level the width of overlapping 
is governed by the mesh size (h and 2h, respectively). Similar variants can be 
constructed for a higher overlapping of subdomains. For the decomposition 
presented in Fig. |4] we have 

i?(Xa) = [0,1/3,2/3,1], a = 1,2,..., p. (77) 

Obviously, in this case we have a grater volume of data exchange, but at the 
same time the transition from one domain to another is much smoother. The 
latter allows us to expect a higher accuracy of the approximate solution. 

Consider now the results of approximate solving problem (j65p - (|67l) . which 
has the exact solution ([74]) . Let ni = 2, n2 — I, T ~ 0.01 and the grid is square 
A^i = iV2. Calculations have been performed using regularized schemes with 
the additive (scheme (O, (gD), (gS]), (gS])) and muhiplicative (scheme 0, (|iT|). 

(|49pl perturbation of the transition operator at cr = 1, as well as vector 
additive scheme ([55|) . The results are compared with the difference 
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solution obtained via implicit scheme dH), ([B]), © at cr = 1. The error of the 
approximate solution was estimated through value e(i") = ||2/"(x) — w(x, 
at a particular time step. 



implicit 
additive 
muitipiicative 
vector 




Fig. 5 Error at Ni = N2 = 32 and N = 10 



Considering decomposition ([75]) (the width of overlapping is h) with the 
space grid Ni — N2 ^ 32 and time grid = 10 (r = 0.001), we can compare 
the error norms of the difference solution obtained using different schemes (see 
Fig. [5]). Figure [BHl] shows the local error at the final time moment. The error is 
localized in the area of overlapping and it is much lower for the vector scheme 
of decomposition in compare with the additive and multiplicative variants of 
regularized additive schemes. 

In contrast to the implicit scheme, the error of the approximate solu- 
tion grows with increasing the space grid for domain decomposition schemes 
(Fig. [Q]). In this case the width of overlapping is reduced by half. 

The dependence of results on the width of overlapping is shown in Fig.[TOl It 
is easy to see that decomposition (j77p demonstrates the approximate solution 
of essentially higher accuracy in compare with decomposition (l75l) (compare 
Fig. E] with Fig.Iini). 



7 Conclusions 



1 . In this paper we have constructed the operators of domain decomposition 
for solving evolutionary problems. Splitting of the general not self-adjoint 
nonnegative finite-dimensional operator is performed separately for its self- 
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Fig. 6 Error of scheme ([Tjl, ||4TJ, igSj, ll49l l 



l.Oi ^ , ^ ^ 1 I — ^0.15 




Fig. 7 Error of scheme (0, Bill . JiSl l. (ggj 



adjoint and skew-symmetric parts. This preserves the property of nonneg- 
ativity for operator terms associated with individual subdomains. 
2. There are constructed unconditionahy stable regularized additive schemes 
for the Cauchy problem for evolutionary equations of first order based on 
splitting of the problem operator into the sum of not self-adjoint nonnega- 
tive operators. Regularization is based on the principle of regularization for 
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Fig. 8 Error of scheme (|33]l, (ESj, ll56ll 
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Fig. 9 Error at Ni = N2 = 64 and N = 10 



operator-difference schemes with perturbation of the transition operator for 
the expUcit scheme. Both additive and multiphcative sphttings are consid- 
ered. It was highhghted the relationship of such regularized schemes with 
the additive schemes of summarized approximation: additively-averaged 
schemes as well as standard component-wise splitting ones. 
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Fig. 10 The error for Ni = N2 = 32, Af = 10 and decomposition R = [0, 1/3, 2/3, 1] 



3. Vector additive schemes of full approximation are selected among the split- 
ting schemes for evolutionary equations. They arc based on the transition 
to a system of similar problems with a special component- wise organization 
for searching the approximate solution at the new time level. 

4. The numerical solution of the boundary value problem for the parabolic 
equation in a rectangle was conducted. The calculations allow to compare 
various schemes of domain decomposition and to show the accuracy depen- 
dence of the approximate solution on the width of overlapping. The vector 
additive scheme of domain decomposition demonstrates the best results in 
terms of accuracy. 
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